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Uniqueness of effective interaction defined in an extension of the Kohn-Sham theory is proved, if 
the model with a non-degenerate ground state exists and to reproduce a correlation function as well 
as the single-particle density of an electron system. The two-body interaction term is regarded as a 
Hubbard-type short range interaction term. The interaction strength is a functional of an element 
of a two-body reduced density matrix of the electron system. The uniqueness theorem gives a basic 
principle for effective description of electron systems. 
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Progress in theoretical description of many-electron 
systems and the high-performance-computational tech- 
niques have allowed us to predict new functionality of 
materials and even to design unknown materials by theo- 
retical model calculations. The density- functional theory 
(DFT)0,HH|1|| has played a central role to realize reli- 
able and powerful skills for the material scientists. More- 
over, DFT is valuable for pure physics since it provides a 
way to create an effective Hamiltonian of an interacting 
quantum system. 

By introducing an idea that the single-particle den- 
sity is reproduced by an effective many-body system, 
Kusakabe derived a kind of extension of the Kohn-Sham 
theory. In our extended Kohn-Sham scheme, we have 
an effective many-body Hamiltonian expressed in second 
quantization. The creation and annihilation operators 
for artificial particles arc defined using single-particle 
wave functions which may be expanded in Kohn-Sham 
orbitals. The effective potential for our modified Kohn- 
Sham single-particle equation is defined with a resid- 
ual exchange-correlation energy functional E IXC [n] . Here 
n(r) is the electron density. This residual term is deter- 
mined via selection of the interaction terms introduced 
explicitly in the model. If an interaction term is written 
in terms of a localized orbital, the interaction may be 
interpreted as the Hubbard-type short-range interaction. 
Shift in interaction strength may change -E^c [n] , but the 
formulation holds without uncertainty. This arbitrari- 
ness means that our effective theory is not uniquely de- 
termined within the framework of the density-functional 
theory. 

This problem of uncertainty is solved in this paper 
by formulating the theory as a density-matrix-functional 
theory (DMFT). DMFT is originally written as an ex- 
tension of DFT. In the usual DMFT, m-body reduced 
density matrices JmM, ••• , r' m ; ri, • ■ • , r m ) are used as 
the basic variable. 0, 0] Here is a space coordinate for 
electrons. Since a two-body reduced density matrix gives 
the exact interaction energy, knowledge of the reduced 
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density matrices determines the true ground state en- 
ergy of a many-electron system. However, there remains 
the 7V-representability problem in the known DMFT. 

On the other hand, a natural form of two-body reduced 
density matrices appears in our extended Kohn-Sham 
scheme. The form representing density fluctuation can 
be a density-density correlation function in a localized 
orbital. We will show that the simplest form is positive 
real- valued and its range is from to 1. Our model gives 
a minimization process of a wave- function functional and 
the minimum of the model is shown to exist. If we re- 
quest that the fluctuation in the original Coulomb system 
is reproduced in our model, which should have a non- 
degenerate ground state, the model of the many-electron 
system is uniquely determined. In a specified situation, 
we will conclude an existence condition of the properly 
defined extended Kohn-Sham model with an interaction 
parameter U as a density-matrix functional. 

We may choose relevant orbitals and introduce an op- 
erator detecting fluctuation per each orbital. To simplify 
the discussion, we consider at first an impurity problem 
with a single relevant orbital (f>i(r) at an atomic center 
in a bulk. This situation naturally gives a Kondo prob- 
lem. However, it is only technical to consider an impurity 
problem. Once our strategy is given for this typical case, 
further possible extension becomes trivial. 

We introduce an operator defining a correlation func- 
tion (n 2 ), which is the density-density correlation func- 
tion on 4>i(r). 

ill 2 ,) = ((n i)T + m,i - n lA - n lA f) (1) 

Here, <pi(r) is normalized and defines associated creation 
and annihilation operators denoted by c\ a and Cj jCr . The 
number operator ni_ a = c| CT Ci j(T is used. Expectation 
values hi_ G = (^\ni_ a \^!} are taken for a many-body state 
| , F), which is a trial state unless explicitly specified. 

The density n(r) and the expectation value (n 2 ) given 
by are the fundamental variables in our theory. For 
the ground state \$g) °f the Coulomb system, ri(?(r), 
(h 2 )g given by |$g) w iU be used. We are seeking for a 
model in which ng(r) and (n 2 )c ar e reproduced exactly. 

We first show the next lemma. 
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Lemma 1 

(n?) is real. The next inequality holds. 

< (n 2 ) < 1 . 



(2) 



Proof: Since rij i(T is Hermitian, is real and positive- 
scmi-dcfinitc. Arbitrary |^) is expanded as, 

|*) = A \y )+Aicl t \y 1 )+A 2 clj* 2 )+A 3 c\ A cl l \*3) . 

Here each (i = 0,1,2,3) satisfies Cj )(7 |*j) = 0. Co- 
efficients Ai satisfies a normalization condition \Ao\ 2 + 
|^i| 2 + \A 2 \ 2 + \A 3 \ 2 = 1. Using Ai, we have 



2|^ 3 | 2 + (|A 1 | 2 



W 2 | 2 + 2|A 3 | 2 ) 



x {l-(\ Al \ 2 + \A 2 \ 2 + 2\A 3 \ 2 )} 



(3) 



This expression has the maximum value 1, only if \Aq\ = 
\A 3 \ = l/y/2 and \A X \ = \A 2 \ = 0. 

We require (n 2 ) to be reproduced in our model system. 
There are some possible manners to formulate the density 
matrix functional theory. One might consider generaliza- 
tion of constrained search of Levy 9] in order to make a 
density matrix functional theory. However, this is not the 
only way. We utilize Levy's universal energy functional 
F[n] given by 

F[n] = mm(*'|f + V„\V) , 

to formulate our theory. Using the electron field opera- 
tors ip^ (r) and ip a (r) , the kinetic energy operator T and 
the Coulomb interaction V ee are given by, 



U no = - ld 3 rd 3 r' 



-^^(r)^(r')^'(r')^(r) 



It is known that a minimizing in the definition of F[n] 
exists, although the convexity of F[n] is denied. |f| 

Now we consider the reduced interaction term, which 
may be regarded as correlation or fluctuation given by, 



(4) 



We define a functional G[$], which is an extension of 
Hadjisavvas-Theophilou's functional. [fil llflj 

G[tf] = <*|T + U rcd |*) - min + V Ted \V') 



+ F[n 9 ] + J d 3 rv oxt (r)n^(r) 

= mf + v Ied \v 



1 r %(r)^(r') . 3J 
2./ |r-r'' 



[n*] + / d 3 rw oxt (r)n*(r) 



(5) 



Here, the residual-exchange-correlation functional is, 

= min (*'|f + y cc |*') - min (*'|f + U rod |*') 
1 /■ ^(r)^(rQ d 3 rd 3 r / 



In general, the definition of V re d determines E rxc . Vary- 
ing U may change E rxc . Any U is allowed to determine a 
model. However, almost all of the models are not proper, 
since fluctuation may not be reproduced. To make U de- 
pendence explicit, we utilize a symbol G[/[W] below. 
We next prove a set of exact statements in our theory. 



Lemma 2 Assume that the ground state of a Coulomb 
system given by v ex t(r) is non-degenerate, i) The ground 
state ^ of a corresponding extended- Kohn- Sham model 
Gu[^\ with a given positive U exists, ii) For fixed n(r), 
F(U) — min*_ m (\f , |T'-|- ^n 2 \ty) is a continuous function 
of U. iii) If a state |^) is the ground state of G[/ 1 [\f r ] 
and Gu 2 [^} with < U% < U 2 simultaneously, \H>) is the 
ground state of Gufif] in a finite range [Ui,U 2 ] ofU. 

Proof: The ground state, which gives the minimum of 
G(j[^], satisfies two conditions. (1) |\&) — > nc(r). (2) 
|*) minimizes (*|T + ^n 2 |*). The first statement i) is 
proved followingthe same argument in the proof of the 
theorem 3.3 in |||. The only requirement which has to 
be shown is that (T + ¥rif ) is a positive quadratic form, 
which is trivial from the definition. 

The second statement ii) is now discussed. Assume 
that F(U) is not continuous when U — Uq > 0. This is 
equivalent to a statement that 3 e > 0, v <5 > 0, 3 C/ > 
0, s.t. {\U -U \ < S & \F(U) - F(U )\ > e}. For 
simplicity, let us further assume that < Z7 — Uo < S 
and F(U) — F(Uq) > e. If we let 8 — e, we have, 

U 



min (M/'|T+^n 2 |vI/') 
> min (*'|f + -y^i I*') + £ 



3" 



= <*o|T+^n 2 |*o)+ £ 
= (* |f+^n 2 |* ) + ^_^ 



*oK 2 |*o) 



u 



s 



> {*o\T+~rg\9 )--+e 
= (*o|T+|n 2 |* ) + |. 

Here we have utilized the lemma ^ This inequality con- 
tradicts the definition of min*/^ n (^'|T + -jn 2 |*'). In 
other possible cases, we have the same conclusion. 

The third statement iii) is shown as follows. We define 
U G = n G^)}- Assume that a |*) mini- 

mizes both G Vl [*] and G U2 [*]. Then for any |*') e Kg 
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we have the next inequality with two parameters a > 
and /3 > 0, which do not simultaneously become zero. 



.^qf + ^n 2 |*') 



> a(*|T- 



Ui 2 



nf|*) + /3<*|T- 



^n?|*> 



This inequality ensures that for {/' = U ^a+pf 



U' 



-n?l*) 



A minimizing |*') of G^[$] is given by |*). Thus I*) 
has to be always a ground state of Gjj [*] , if t/ satisfies 
U 1 <U<U 2 - 

We now define a concept of "a proper model" . In the 
present context, a model is properly determined, if its 
lowest level is non-degenerate and if it reproduces nc(r) 
and (nf ) G at the same time. The second condition should 
be changed, if the definition of Vr e d is changed. The cen- 
tral problem of the present study is summarized in the 
next question. If we find a proper model Gj/f*] by op- 
timizing U, is there another U' which do give a proper 
model? Our uniqueness theorem given below denies ex- 
istence of the other model. 

Suppose we found a proper model with U = U\. The 
ground state |*) of Gui [*] is non-degenerate and thus 
unique. If we shift U, |*) may found as a ground state of 
G;j 3 [*] with e.g. U 2 > U\, in which the lowest level may 
be degenerate. Lemma |2 hi) tells us that, in a finite re- 
gion [Ui, U 2 ], |*) is always the ground state of the model 
£rj/[*]. The ground state |*) of the model does not feel 
the short range interaction by U. We call the situation 
"the irrelevant U case" This happens when fa is out of the 
system with no electron occupying fa or fa is occupied 
by a single electron without hybridization. Another case 
is when the system shows the complete ferromagnetism. 
We do not intend to consider these situations. In real 
materials, we have finite hybridization matrix elements 
between fa and a neighboring orbital. Without the large 
magnetic field, the true complete ferromagnetism never 
exists in real solids. Thus we assume that |*) does not 
appear as a ground state of G(y[*] except for U = U\. 
We do not assume, however, that G(y[*] may have de- 
generacy at U ^U\. We now state our main theorem. 

Theorem 3 Assume that the ground state of a Coulomb 
system is non- degenerate. A proper extended- Kohn- Sham 
model given by Gjj[^] which has a non-degenerate ground 
state and reproduces both n G (r) and (nf) G is uniquely 
determined, or it does not exist. 

Proof: Assume that we have U\^U 2 and two extended 
Kohn-Sham models, i.e. G^ [*] and Gj/ 2 [*], reproduce 
na(r) and (t^)g- Because we do not consider an irrel- 
evant U case, l^Ui) an d |*c/ 2 ) are different with each 
other. Let us define E^ s be the ground-state energy and 



Erxc[na] be the residual exchange-correlation energy of 
these models, (i = 1, 2.) We have the next inequality. 



E, 



(2) 
GS 



(*u 2 \f+^n?\* U2 )+Egl[n G ] 
e " fji ^3„/^G(r)n G (r') f ~ 



■ > J ^dV '^^j' ' + j tfVu cxt (r)n G (r) 
(*u 2 \T + ?±rg\9 Ua ) + £«[n G ] 

\r-r'\ J 

+ \(U 2 - UxWut \n*\* U2 ) + E£l[n G ] - £«[n G ] 

^,3, n G{r)n G W) f , 3 



— / tfVdV 
2 ./ v V 



+ / cPrv CKt (r)n G (r) 



\{U 2 - UJiVui In?!*!,,) + E£l[n G ] - Eg\[n G ] 

<*aJT + ^n 2 |* ai ) + i4 2 >[n G ] 

e " ^3^G(r)n G (r') , f , 3 



, cfVtfV 
2 7 r - r' 



+ / d d rv cxt (r)n G (r) 



This contradicts that |*j, 2 ) is the ground state of a 
model Gjj 2 [*] ■ The above argument does not deny non- 
existence of such a model. 

Thus, we have shown that U is uniquely determined, if 
a proper model exists. Does U really exist? If we have a 
case (called a case II below) in which Gjj [*] has always 
a unique ground state irrespective of U, we can show a 
sufficient condition for existence of the proper model. We 
have a next corollary. 

Theorem 4 Assume that the ground state of a Coulomb 
system is non- degenerate and that we have a case II. 
An expectation value (n 2 ) given by the minimizing |*) 
of Gj,[*] is a function of U. The map from [0, 00] to 
[0, 1] is one-to-one, continuous and thus monotone. 

Proof: Since we have Lemma^ we have to show that the 
map is a one-to-one mapping at first. Let us assume that 
we have two J7,s (i = 1, 2), satisfying < U\ < U 2 < 00, 
for which the minimizing |*i/ 4 ) of Gj/J*] reproduces 
n G (r) and expectation values (n 2 ) at the same time. Us- 
ing the same inequality used in the proof of Theorem |3| 
we see that the first assumption is in contradiction to the 
assumption of the case II. 

On the continuity, let us assume the contrary. Then 
we have at least one Uq on the U axis, where the discon- 
tinuity occurs. In another word, 

3 e > 0, 3 [/. i > (i = 1, 2, • • • 00), s.t. { lim U t = U 

i—^oo 

&|(*!/ i |n 2 |*[, i )-(*c; |n 2 |*t, D )|>£}. 
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This condition means that \^Ui) 7^ l^f/o)- Now, the con- 
tinuity of F(U) shown by Lemma 2 ii) and the selection 
axiom ensures that 

v e>0, \ s.t.\F(Ui) - F(U )\ <e. 

Namely, \^Ui) is also the ground state of G(Uq). How- 
ever, since \^Ui) \^u ) f° r au ^ii this leads to a con- 
tradiction to the uniqueness of \^u ) assumed in the case 
II. Immediately, monotonous behavior is concluded. 

Now we define a set of numbers R ( p i = {r\r € 
[0, l], 3 U s.t. r = (^ulllil^u)}- We have immediately a 
next corollary due to Theorem 0] 

Corollary 5 There exists a one-to-one onto-mapping 
from R t f )i onto [0, 00], which is defined by an inverse map 
of the function of U in Theorem ^[ 

We define a way of modelling of the Coulomb system 
as a process to find G;y[ , I'] by determining fa and U. We 
have immediately a next physical theorem. 

Corollary 6 // and only if (n^c * s in R<t>a a proper 
extended Kohn-Sham model exists. 

Proof: The statement is trivial due to Theorem [S] 

Here wc define three minimizing states, I'i'o), V&u) and 
|*c), given by 

min = (* |T|*o) , 

Wo — mq, 

min = (*u\f + -Tg\*u) , 
min = (9 c \f + V ee \*c) ■ 

Apparently, ^ c = * GS . Wc have (nf) = (*o|^|*o>, 
(n?)u = (*c/|nf|*c/), (nf)c = (*cs|nf |*gs>< These 
definitions and statements above conclude the followings. 

Lemma 7 The following holds. 

1- (™?)e/=oo = min*^ na ,(n?). 

2- (nf)os > (ni)u=oo. 

Due to Theorem 0| we can say the following condition 
for U to exist. 

Theorem 8 Consider a case II. Then, R r j >i — 
[(Hi)(7=oo, (nf)o]- If (n%)c < (n?)o, the proper extended 
Kohn-Sham model exists and the unique U given by The- 
orem 3 exists. 

Now we discuss several remaining problems. We may 
search for fa to find a proper model. However, we should 
note that both (tj|)g and R ( p i shift, if fa is changed. For 
a lattice, we may consider 

(V red > = §l>* 2 >- (6) 



The same argument as the above is possible. Then one 
would find U[(n*) G ]. 

To determine fa, we may see the value of (n^a- Prac- 
tically, we can check orbital energy Ei of fa , which is an 
expectation value of the one-body part of the extended 
Kohn-Sham equations. @ We have a Fermi level Ep ten- 
tatively at the step to obtain the Kohn-Sham orbitals. If 
£i is far below Ep or far above Ep, the occupation of <pi 
becomes 2 or 0. Then, we inevitably have (nf )g = 0. If e, 
is around Ep and if the system is in a delocalized scheme, 
= |4| = |4| = |4| - 1/2 and (n?) G ~ 1/2. How- 
ever, if the Coulomb correlation make reduction in (nf )g 
for £i is around Ep, we are required to introduce a li- 
teral to reproduce the density-density correlation func- 
tion. This guide line can be a determination rule of min- 
imal numbers of localized orbitals fa, for which we need 
to set [/-terms up. If singly occupied configuration may 
not appear in the true ground state, then, |4| = 42 1 — 
and (nf)G — 1. This condition may be given if the state 
is superconducting. In this situation, negative U may be 
required, although we cannot utilize Lemma [21 i) for this 
case. 

Lemma |5] iii) leads us to four possible scenarios for 
G[/[ V I'] in the whole range of < U < 00. Here we do 
not consider the trivial degeneracy coming from SU(2) 
symmetry of the total spin and take the highest weight 
state as a representative state. 

Case I A single state \^>) is always the ground state of 

Gum. 

Case II For each U, a unique state |\&) different from 
others appears as the ground state of Gjj [*£] . 

Case III In a finite region [U\, U2] of U, a state is a 
ground state of Gy^]. For any other U, the state 
1^) never becomes the ground state. 

Case IV For any Uq, a ground state of Gu [^] never 
becomes a ground state of G^^] with U 7^ Uq. 
There is a U , for which two (or more) states become 
the ground state of G\j [&] . 

The case I is an irrelevant U case. If \^f>) reproduces (nf) 
in a case III, this is an irrelevant U case. In the case III 
or in the case IV, we have at least one phase transition 
with respect to U in the model. If {t^)g is reproduced 
by a ground state of a degenerate model in a case IV, we 
have to deny the model, since the degeneracy is artificial 
and the model is improper. 

If we do not have a case II for a Coulomb system, 
although F(U) is continuous, (rif) given by \^u) may not 
be unique because of possible degeneracy and may have a 
jump as a function of U. The value of {nf)c might be in 
this jump. In this case, the model has an artificial phase 
transition and thus the model is not proper. Therefor, 
we have to search for a proper model with transferability 
against change in environment by changing fa. If there 
exist discontinuity in (n?) along the U axis, it is detected 
by a jump in (T), which is concluded by Lemma JSjii). 
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Finally, we would like to comment on conceptual 
framework of the whole study. You may think that the 
case II may be special. Arbitrariness remains in the se- 
lection of (pi. However, we can say that Theorems|3]and[S] 
are new guiding principles to construct an effective the- 
ory of complicated correlated systems. In our model, 
we may choose only one [/-term, if appropriate. This 
is impossible for a method to make an effective model 
using reduction of the whole phase space by choosing a 
seemingly proper basis. Our method is completely dif- 
ferent from the renormalization group approach. In our 
method, we keep the property that nc(r) is always re- 



produced and obtainable. Otherwise, an inverse problem, 
e.g. the materials design, would not be solved properly. 
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